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Abstract. - We apply the generic dynamical scaling theory (GDST) to the surfaces of CdTe 
polycrystalline films grown in glass substrates. The analysed data were obtained with a stylus 
profiler with an estimated resolution lateral resolution of Ic — 0.3 /im. Both real two-point 
correlation function and power spectrum analyses were done. We found that the GDST applied 
to the surface power spectra foresees faceted morphology in contrast with the self-affine surface 
indicated by the local roughness exponent found via the height-height correlation function. This 
inconsistency is explained in terms of convolution effects resulting from the finite size of the probe 
tip used to scan the surfaces. High resolution AFM images corroborates the predictions of GDST. 
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The kinetic roughening of interfaces is an outstanding 
topic of nonequilibrium Statistical Physics which has been 
intensively investigated in both theoretical [l}|5| and ex- 
perimental 6^10 frontlines. A generic dynamical scaling 
theory (GDST) for evolving interfaces includes both in- 
terface fluctuations and power spectra (structure factor) 
in the real and momentum spaces, respectively The 
analysis in the real space can be performed using the local 
height-height correlation function 

(1) 



G{l,t) = {[h{x + l,t) -h{x,t)]^), 



where the overbar means averaging over the surface profile 
and (• • •) the averaging over distinct profiles. GDST gives 
that G{l,t) = P°'^{l/t^/'') where a and z are the rough- 
ness and dynamical exponents, respectively. The scaling 
function $ behaves as 



-la 



if X > 1 



(2) 



The local roughness exponent aioc determines the scal- 
ing locally. The power spectrum, defined as S{k,t) — 
(/i(k, t)ft,(— k, t)) with h being the Fourier transform of 
the height surface, scales as S{k,t) = k-^'^°'+'^'>'^{kt^^^), 
where the scaling function 5" is 



■9{x) 



,2(c 



if a; > 1 



(3) 
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In this scaling function, d is surface topological dimension, 
and as is the spectral roughness exponent. 

This generic scaling ansatz has been observed in a large 
collection of models and experiments as can be looked 
up in (l][6|[8] and references therein. This scaling ansatz 
implies a constraint between the exponents that depends 
specially on the spectral roughness exponent [l|[2,ll . If 
< 1 the surface is self-afhne with spectral and local 
roughness exponents being equal, aioc = «s- If ckj > 1 the 
surface is locally smooth with aioc — 1- Each case is still 
classified into two subclasses. For as < 1, if = a we 
have the regular Family- Vicsek (FV) scaling, otherwise 
the system has the intrinsically anomalous scaling. For 
ag > I, we have the super-roughening scaling ii a^ = a 
and faceted growth scaling otherwise. The constraint be- 
tween as and aioc is a beautiful analytical result derived 
from the scaling ansatz 



11 while the subclasses are al- 



lusive to their physical implications fH , with exception of 
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the FV scaling which is due to the ansatz conceivers 

FV, super-rough, and intrinsic scalings have been re- 
ported in several experimental works [l][6][8] but the 
faceted one has only been achieved quite recently in the 
electrodissolution of pure polycrystalline iron [6]. The 
surfaces undergo a transition from intrinsic to faceted 
anomalous scalings as dissolution time increases. The 
surfaces in the faceted regime were highly anisotropic 
and, consequently, the scaling analyses were performed 
with one-dimensional profiles, in a direction orthogonal 
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to the anisotropy, in addition to the analyses of the two- 
dimensional surfaces. The faceted anomalous scaling was 
evident only for the one-dimensional case. This phe- 
nomenon was ascribed to the averages over all directions 
performed in d = 2, which underestimate the scaling ex- 
ponents due to the contributions of unclear faceted mor- 
phology in some directions. 

We investigate the anomalous roughening in semicon- 
ductor CdTe polycrystalline thick films grown on glass 
substrates. The real space scaling of this system was previ- 
ously investigated [9|[T3). In the present work, we identify 
a crossover in the dynamical scaling exponents not no- 
ticed in Ref. [9| and perform the power spectrum scaling 
analysis. We find that the GDST applied to the surface 
power spectra foresees faceted growth whereas the correla- 
tion function analysis in points out for a self-affine surface. 
This inconsistency is explained in terms of the convolution 
effects resulting from the finite size of the probe tip used 
to scan the surfaces. 

The CdTe films were deposited on glass substrates cov- 
ered with a transparent conducting oxide (TCO), Sn02:F, 
resulting a rough initial surface. The glass/TCO surface 
has a rms width of 16 nm that is much larger than 3 nm 
observed for the pure glass substrate. Although the pres- 
ence of this layer introduces an undesirable large initial 
roughness, it is important from an application viewpoint 
because it is the front contact layer of solar cells produced 
with this material 
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The samples used in the present work were produced 
by hot wall epitaxy at a growth rate of approximately 0.1 
nm/s. Details of the growth system and sample prepa- 
ration can be found elsewhere [SjjlS]. It is important to 
mention that the growth of some samples was repeated 
using the molecular beam epitaxy technique and they ex- 
hibited the same behaviour. The scaling analysis was done 
in one-dimensional profiles with 300 /im of length and 4570 
pixels acquired with a stylus profiler XPl-AMBIOS. Con- 
sequently, the topological (scan) dimension is d = 1 in- 
dependently of the 2 -f 1-dimensional surface. At least 
20 profiles scanned at random directions were used in the 
averaging. 

Despite of the good vertical sensibility (better than 1 
nm) and step size scan (better than 100 nm) of the profiler 
the large tip radius (« 1 /im) limits the lateral resolution 
power due to convolution efi'ects [16]. As a consequence, 
the data at short scales correspond to the convolution of 
the film surface and the probe tip geometries. A criterion 
accepted by the microscopy community establishes that 
a probe tip can resolve peaks separated up to 20% of its 
radius. Therefore, one can estimate a lateral resolution 
Ic ~ 0.3 jjLxa. for this scan device. This lateral resolution, 
however, does not hinder the measurement of long wave- 
length fluctuations and the global roughness of the entire 
profile can be accurately determined. The short scale de- 
tails of the surface morphology were visualized using an 
atomic force microscope (NTMDT-Ntegra Prima) operat- 
ing in the semi-contact mode. These AFM images were 



not used in the scaling analyses. The use of profiler scans 
allows to study large amplitude and wavelength fluctua- 
tions not achievable with a regular AFM. The growth time 
was varied from i = 30 to 660 min and the substrate tem- 
peratures T = 150, 200, 250, and 300 °C were analysed. 
As will be shown, reliable scaling of the power spectrum 
were obtained for T = 250 and 300 °C. The CdTe films de- 
tach from the glass substrates and crack during the cooling 
process for growth times longer than 660 min. This tech- 
nical limitation hindered the growth of samples beyond 
this limit. 

The exponents aioc and as were obtained directly from 
the scaling behaviours of G{l,t) ^ p'^ioa q^^^ S{k,t) ^ 
j^-{2a,+i) ^ ^ ggQ rpj^g global interface width is 

defined as 



w{t)^([h{x,t)-h{t)Y 



1/2 



(4) 



where /3 is the growth exponent. GDST states aioc = 
min(l, as) and a = /3z p|[l7]. Even though the local and 
spectral exponents are not independent, the local rough- 
ness exponent was measured for a comparison between the 
theoretical prediction and experimental observation. 




Fig. 1: Interface width evolution for two growtli temperatures. 
Lines are power law regressions. Inset shows two surface pro- 
files for growth times t — 180 min (dark) and t — 660 min 
(light) at T = 300 °C. 

The surface width evolution, which is shown in Fig. [l| 
has a crossover from a regime of low to high growth ex- 
ponents. This crossover was not perceived in the former 
scaling analysis of this CdTe/TCO/glass system, but it is 
also present as one can check in Ref. [o]. The crossover 
is probably on account of the initial surface roughness, 
which initially diminishes due to the occupancy of grooves 
and valleys. The temperature increase enhances the down- 
ward funnelling, reducing the crossover times and, conse- 
quently, the largest power-law interval was obtained for 
T — 300 °C. An additional evidence for our proposition is 
that this crossover was not observed in the CdTe deposi- 
tion on pure glass substrates (a smaller initial roughness) 
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at the same experimental conditions 15 . Therefore, the 
growth exponents in Ref. [oj are under-estimated but the 
main conclusion of that work, stating a j3 exponent de- 
pending on temperature, still holds. 

Whole scaling analyses were performed for the growth 
times after the crossovers. The /3 exponents are shown in 
Table [T] The large uncertainties for temperatures 150 and 
200 °C reflect the short time intervals of power law growth. 
Taking into account the uncertainties, these exponents are 
nearly constant for temperatures 150, 200, and 250 °C. For 
T — 300 °C, the growth exponent /3 > 1 is a signature of 
an unstable growth that is unusually fast in the framework 
of kinetic roughening. However, self-affine (FV scaling) 
unstable growth was recently found in stochastic equations 
related to nonlocal interface dynamics (3[|4|. 

Interestingly, the growth exponent obtained numeri- 
cally for the unstable version of the stochastic Michelson- 
Sivashinsky (SMS) equation, Psms = 1-14 [s], coincides 
with our estimate /? = 1.11(2). The number in paren- 
thesis represents the uncertainty in last digit meaning 
j3 = 1.11 ± 0.02. The mentioned unstable SMS equation, 
written in the momentum space with suitable rescaled pa- 
rameters, reads as [sj 



dth = (/c - k^)h + ^J^(|V/i|2) + f] 



(5) 



where J-{g) = g is the Fourier transform of g and "q a 
Gaussian noise. Apart from the term fc/i, this equation is 



the well known Kardar-Parisi-Zhang equation 18 related 
to a non-conserved growth. The instability is due to the 
term kh that introduces the non-locality in the equation. 
However, as shown in the present work, the CdTe/Glass 
system is not self-affine (a^ ^ aioc) and cannot be directly 
associated with the SMS equation. 



Power spectra for T — 300 °C are shown in Fig. 2(a) 



The curves exhibit the usual plateaus for small and power 
law decays for large momenta. The power law yields 
as = 1.27(5) and, consequently, GDST states aioc = 1- 
Notice that momenta in the range k < 2tt/Ic were used 
in this analysis. Moreover, the curves shift upwards for 
large k as time increases, representing a > as since 
S{k,t) ^ ;,2Q,+i^2(a-Q^)/z ^Yie scaling ansatz given 
by Eq. ([s]). GDST thus predicts surfaces with anoma- 
lous scaling and facets. Even though the profiler provides 
coarse-grained scans that hinder to reliably resolve faceted 
morphologies, as can be seen in Fig.jga), a AFM image 
with a scan step size of 5.8 nm and probe radius of 10 nm 
and, therefore, a resolution about 100 times better than 
the profiler, clearly exhibits the faceted surface morphol- 
ogy predicted by GDST, as shown in Figjsjjb). That is 
a remarkable feat of the GDST which has foreseen the 
faceted morphology in coarse-grained data with a poor (if 
any) resolution of the facets. From a theoretical point 
of view, this result corroborates that GDST connects the 
scaling properties of long wavelength fluctuations with lo- 
cal morphological properties. In our case, the scaling prop- 
erties of mesoscopic {2tt/Ic < fc <C 2tt/L) and macroscopic 
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Fig. 2: (a) Power spectra of scanned profiles for T — 300 °C. 
Growth times increase from bottom to top. The dashed line 
has a slope -3.6 and represents the scaling S ~ (b) 
Collapses of the power spectra using exponents a = 1.94 and 
z — 1.75. Dashed lines represent the asymptotic behaviours of 
the scaling function ^{x). The slopes are 2a -I- 1 ~ 4.9 and 
2(a - as) ~ 1.3. 



(fc <C 2tt/Ic) scales are connected with the surface mor- 
phology in a microscopic scale fc > 2tt/Ic- 

In Fi^ 2(b) I we probe the scaling ansatz for "^(x), given 
3), by plotting k'^°'+^S{k,t) versus kt^/'-, a = (iz, 



by Eq 

where z = 1.75 was obtained via gradient correlation 
function [ explained in the next paragraph and shown in 
4(a)] and /3 = 1.11 via global interface width (Fig- 
The excellent collapse obtained for t > 90 min 



Fig. 



provides an additional evidence that the data is described 
by anomalous scaling for faceted surfaces. Notice that the 
curve corresponding to t = 30 min, which was removed 
from the determination of the scaling exponents (Fig. [T]), 
does not collapse. Finally, the asymptotic scaling forms of 
^(x) are also confirmed as indicated by the dashed lines. 

The dynamic exponent z can be determined using 
the slope-slope correlation function T{l,t) = (V/i(a; + 
l)Wh{x)). The correlation length ^, which is assumed to 
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(a) 



(b) 



Fig. 3: Surface morphologies for CdTe after a growth time of 
300 min at a temperature of 300 °C. (a) A 10 ^m scan with 
using a profiler, (b) A 3 /im x 3 /im AFM image illustrating 
the faceted CdTe surface. 



scale as C ~ t^/^ in GDST, can be defined as the first zero 
of the correlation function. The correlation functions and 

There is no clear power 



lengths are shown in Fig. 4(a 



law regime with f > for T === 150 - 250 °C. For 300 °C, 
the slope provides \/z — 0.57. 

Power spectrum analysis is itself sufficient to carry out 
the GDST, but correlation functions are widely used in 
the experimental investigations. Fig. |4(b)] shows that the 
height-height correlation function, equation Q, also fol- 
lows the usual qualitative behaviour of anomalous scaling, 
in which the curves are shifted upwardly as time evolves. 
However, a power law regression in the linear interval in- 
dicated turns out a self-affine surface with a/oc = 0.82(2), 
in disagreement with GDST. The scaling exponents for 
the other temperatures are shown in Table [T] The under- 
estimated values of aioc are reflecting the coarse-grained 
resolution of the scanning device since, in this scale, the 
analysed data are the convolution of surface and probe tip 
morphologies. The inset of Fig. |4(b) [ shows the correlation 
function collapses using the same exponents that collapsed 
the power spectra, as shown in Fig. |2(b)[ The collapses 
are so or more convincing than those obtained for power 
spectra. 

The power spectra for T = 250 °C are shown in 
Fig. 5(a) The curve for t = 600 min also exhibits a 
power decay in a range within k < lix jl^ that yields 
as = 1.6(1). Again, S{k,t) shifts upwardly as time in- 
creases implying a > as- Gradient-gradient correlation 
function does not provide a reliable estimate of the the 



dynamical exponent z [Fig. 4(a) and thus a criterion of 
best collapse for t > 90 min was used instead to obtain 
z = 3.4(4) and a — 2.2(2). The collapse using these expo- 
nents are shown in Fig. |5(b)[ The power spectrum analy- 
ses for T = 150 and 200 °C did not provide reliable scaling 
regimes within the interval k < 2tt/Ic. Even though the 
scaling for T = 250 °C is less precise than for T = 300 °C, 
it still points out a faceted anomalous scaling which was 
also confirmed in AFM images. 

The power spectra for T — 150 and 200 °C do not ex- 
hibit scaling properties accurate enough to determine the 
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Fig. 4: (a) Slope-slope correlation functions for the CdTe sur- 
faces. In the main plot, the correlation for distinct times and 
T — 300 °C are shown. Inset shows the correlation length 
against time. Straight line is a power law regression, (b) 
Height-height correlation functions for T = 300° C. The in- 
terval of regression is indicated. Inset shows the collapse of the 
correlation functions using a = 1.94 and z — 1.75 



exponents but the qualitative behaviour of faceted anoma- 
lous scaling, > 1 and a > Og, were also verified. AFM 
images again confirms the faceted morphologies for these 
temperatures. Indeed, the effect of temperature is to in- 
crease the size of faceted grains. For this reason, the scal- 
ing analysis with profiler date is more accurate for higher 
temperatures. 

Now, we compare the the GDST in the momentum 
space with the results for real space of this CdTe-TCO- 
glass system reported in Ref. 



13 . The scaling exponents 



a and z shown in Table [T] are in disagreement with those 
in Ref. 13 , in which two sources of errors are present 
in the exponent determination. The first one is the un- 
derestimation of the growth exponent used to determine 
a via scaling relation a = /3z. The second one is the 
correlation length ^ that overestimated the dynamical ex- 
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— 30 min 

— 180 min 

— 300miri 

— 660 min 



Table 1: Scaling exponents for different growth temperatures. 
The numbers in parenthesis represent uncertainties in the last 
digit. The aioc exponents were obtained using the height- 
height correlation function given by Eq. The missing 
exponents are due to the absence of a reliable scaling regimes. 



T (°C) 


/3 


as 


z 




150 


0.59(9) 






0.73(2) 


200 


0.56(9) 






0.76(3) 


250 


0.65(2) 


1-6(1) 


3.4(4) 


0.77(2) 


300 


1-11(2) 


L27(5) 


L75(5) 


0.82(2) 



(a) 




(b) 



Fig. 5: (a) Power spectra of scanned profiles for T = 250 °C. 
Growth times increase from bottom to top. The dashed line 
has a slope -4.1 and represents the scaling S ~ fc"'^"'"*"^'. (b) 
Collapses of the power spectra using exponents a = 2.2 and 
z = 3.4. 



ponent z, as we clarify in this paragraph. The correlation 
length was defined in Ref. 



13 as the characteristic length 



of the decay in a two-point correlation function r(/,t) = 
PT[\h{x + l) — h(x)\ < to], where m = 0.1\hmax — hmin\ and 
Pr[A] is the probability that the condition A is satisfied. 
The usual correlation function F = {h{x + l)h[x)) under- 
goes exactly the same effects. The correlation length was 
finally obtained by solving Jq Pdl = 0.1 J^Tdl, where 
T(l) is a two-exponential fit to the experimental data. 
However, if the tail corresponding to the long wavelength 
height fluctuations is or not left out of the regression, the 
exponents may change considerably. In Fig. [6j we compare 
the regressions of the correlation functions in two ranges: 
discarding (range 1) and including (range 2) the tail. The 
regression in the range 2, the same used in Ref. [l3], mis- 
fits the data only for I < 0.3 jim whereas the regression 
in the range 1 fits very well the small scales but devi- 
ates for / > 2 /im, which was the reason for the choice 
of range 2 in Ref. [l^ . The insertion to Fig. |6] shows the 



correlation lengths against time and the respective power 
law regressions. Both cases yield quite satisfactory scaling 
laws, ^1 ~ and ^2 ^ corresponding to ranges 1 
and 2, respectively. Assuming ~ i^/^', where i — 1,2, 
we found zi = 1.72 and Z2 — 3.69. Notice that the re- 
gressions provided two characteristic lengths of the same 
magnitude (1 < £,2/£,i ^ 2) but exhibiting very different 
scalings with time. Repeating this analysis for lower tem- 
peratures (T = 150 - 250 °C) we found < 0.5 /xm and 
no clear power law could be identified as observed for the 
gradient correlation function. Therefore, range 1 is the fit 
interval that yields dynamical exponents consistent with 
GDST. 




l(Hm) 



Fig. 6; Height-height correlation function of the CdTe surfaces. 
The growth time and temperature were 90 min and 300°C. Cir- 
cles represent the experimental data while solid lines the two- 
exponential regressions in the indicated intervals. Inset shows 
the correlation length obtained for the different regression in- 
tervals. 



The conclusions of Ref. 13 are thus partially modi- 



fied since we cannot apply GDST with those exponents 
to infer about which dynamical scaling regime the system 
belongs to 



13 



was 



The exponent z presented in Ref. 
misleadingly associate to the the dynamical exponent of 
the GDST. It would be better to refer to this exponent as 
a long wavelength coarsening exponent, which increases 
with temperature whereas the actual dynamical exponent 
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decreases. These results also stress out the difficulties usu- 
ally observed in the fitting procedure of experimental data 
and the problems which can arise from a misinterpretation 
of the parameters obtained. In this scenario, GDST using 
momentum space appears as a robust method since local 
scaling properties can be obtained from mesoscopic and 
macroscopic scales. 

In conclusion, we presented the scaling analysis of CdTe 
poly crystalline surfaces grown in glass substrates using a 
generic dynamical scaling theory (GDST) [2]. We investi- 
gated both height-height correlation functions and surface 
power spectra. The surfaces were scanned with a stylus 
profiler a resolution estimated as Ic ~ 0.3 /im. 

Surface power spectra exhibit reliable scaling properties 
for I > If. &t temperatures 250 and 300 °C, with a spec- 
tral roughness exponent as greater than 1 but different 
from the global roughness exponent a. The GDST foresees 
therefore faceted morphologies with local roughness expo- 
nent aioc = 1 that represents locally smooth surfaces. The 
height-height correlation function exhibits scaling only for 
I < Ic and the local roughness exponents are in the inter- 
val 0.7-0.8, indicating self-affine surfaces as described in 
a previous study of this system 9 . However, high reso- 
lution AFM images confirm that the surfaces are faceted. 
The underestimated values of aioc are reflecting the low 
resolution of the scanning device since, in this range, the 
analyzed profile is a convolution of the surface and probe 
tip. 

It is worth to note that in a previous reports on anoma- 
lous scaling in faceted morphologies, investigated in the 
dissolution of poly crystalline iron 6! , the faceted anoma- 
lous scaling was evident only for one-dimensional sections 
in a direction orthogonal to the anisotropy. Differently, 
our scaling analysis of polycrystalline CdTe consists in av- 
eraging over several randomly selected directions. 

A renormalization group analysis of stochastic equations 
shows that the anomalous scaling cannot be present in lo- 
cal growth models [l]. Therefore, this theoretical result 
implies that disorder and/or nonlocal effects are respon- 
sible by the anomalous scaling in experimental systems. 
Indeed, the CdTe/TCO/Glass system undergoes both ef- 
fects simultaneously. The amorphous glass substrate re- 
sults in random growth orientations and, consequently, a 
polycrystalline film is obtained. Moreover, the glass sub- 
strate and CdTe films have distinct coefficients of thermal 
expansion implying strained film. The interplay between 
nonlocal strain and substrate disorder is a possible mech- 
anism involved in the anomalous scaling. 

Our work opens new paths to investigate anomalous 
scaling experimentally. A natural extension of this work, 
is to check the effect of the substrate on the scaling prop- 
erties. Scaling analysis of AFM images may also com- 
plement the present investigation, in particular at lower 
temperatures where the scaling properties did not allow 
to determine the scaling accurately. 
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